Identifying contact effects in electronic conduction through buckyballs on silicon 
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We present a theory of current conduction through buckyball (Ceo) molecules on silicon by 
coupling a density functional treatment of the molecular levels embedded in silicon with a non- 
equilibrium Green's function (NEGF) treatment of quantum transport. Several experimental vari- 
ations in conductance- volt age (G-V) characteristics are quantitatively accounted for by varying the 
detailed molecule-silicon bonding geometries. We identify how variations in contact surface mi- 
crostructure influence the number, positions and shapes of the conductance peaks, while varying 
separations of the scanning probe from the molecules influence their peak amplitudes. 
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Molecular electronics represents an ultimate dream for 
nanoscale material and device engineering. Along with 
well characterized, reproducible experiments, quantita- 
tive models for molecular conduction are crucial for a 
proper understanding and benchmarking of this emerging 
field, and for the exploration of novel device paradigms. 
A persistent problem has been incomplete knowledge of 
metal contact microstructures and their influence on con- 
duction. In this respect, a semiconducting substrate pro- 
vides a superior test-bed due to its well studied surface 
chemistry for transport measurements 0,0. It is there- 
fore worthwhile to develop and refine our knowledge base 
using a familiar molecule bonded on a well-characterized 
silicon substrate that leaves little wiggle room for theory. 

Among various molecules probed using scanning tun- 
neling spectroscopy (STS), buckyball (Cqo) molecules 
stand out for their unique well calibrated bandstruc- 
ture, alkali metal doped superconductivity, switching and 
optoelectronic properties Although STS studies of 
buckyballs on metals have allowed detailed comparison 
with theory they do not reveal much information 
about their underlying contact microstructure. In con- 
trast, buckyballs on silicon exhibit considerable varia- 
tion in their G-V characteristics depending on the na- 
ture of their covalent bondings with the surface dimers 

In this paper, we explore conduction through bucky- 
balls on silicon, and correlate observed variations in their 
G-Vs with variations in their contact bonding geome- 
tries (Fig. HJ. A variation in the nature of the molecule- 
substrate bonding leads to a variation in the number and 
shapes of conductance peaks, while a variation in the tip- 
sample tunneling gap leads to a variation in the relative 
peak heights. Our theoretical formulation thus serves 
a dual purpose: on one hand, it tests our quantitative 
model for molecular conduction, in particular on a so- 
phisticated semiconducting substrate; on the other hand, 
it provides useful insights that allow us to deconstruct the 
role of contact geometry on molecular conduction. 

Theoretical technique. We calculate molecular conduc- 
tion by coupling an electronic structure calculation for 
the molecule and the contacts with a treatment of quan- 
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FIG. 1: Different STS measurements 0,E3 on ^60 molecules 
docked onto Si(100) 2x1 surface (left panel). We attribute 
different bonding geometries (middle panel) for each experi- 
ment, leading to a theoretical G-V (right panel) that agrees 
quite well with the corresponding measurement. The arrow 
represents the dimer direction of the reconstructed surface 
(cross going into the page) . The upper geometry corresponds 
to Cqo physisorbed on four surface dimers, the middle one 
represents the buckyball chemically bonded with a single sur- 
face dimer, while the bottom one has the molecule lowered 
into the trough caused by a missing dimer. 

turn transport using NEGF 0, ^} . The Ceo structure 
and Hamiltonian are obtained using density functional 
theory within the local density approximation (LDA). 
The reconstructed surface geometry of Si(100) is ob- 
tained by LDA optimization with a norm-conserving 
pseudopotential in a plane- wave basis [ljj. Although 
it is possible in principle to describe the silicon band- 
structure using DFT, as has been customarily done in 
the past for metal substrates [HI , adapting the same 
process to semiconductors is quite challenging, given the 
complicated bandstructure, extended band-bending and 
incomplete screening, reconstruction, and surface chem- 
istry of silicon. Fortunately, within the NEGF formalism 
one can formally partition the problem so that the only 
quantum effect of the silicon substrate that the molecule 
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FIG. 2: Calculated surface bandstructure of asymmetric 
dimer (solid lines) and symmetric dimer (dashed lines) recon- 
structions along Si(100)-2xl using EHT bulk Si parameters 
[lH and optimized surface geometries [TH. 

is sensitive to resides in its surface Green's function. In 
the past, we developed a technique for combining differ- 
ent electronic structure codes by matching their interfa- 
cial Green's functions expressed in two different atomistic 
basis sets. The match is exact for two equally complete 
basis sets and a best-fit otherwise, and only assumes local 
separability of their one-electron potentials 0, • 

We use an Extended Huckel (EHT) type model 
parametrized by Cerda et. al. [3 to generate a good 
quantitative description of the bulk silicon bandstruc- 
ture (calibrated with LDA+GGA calculations), and also 
the surface band structure of 2x1 reconstructed Si- (100) 
(calibrated with experiments and GGA calculations |l5j). 
To check the properties of Si(100)-2xl reconstructed sur- 
face properly, a slab of 13 silicon layers with a hydrogen- 
passivated bottom layer is used to simulate the band- 
structure of Si(100)-2xl asymmetric dimer (AD) recon- 
struction and symmetric dimer (SD) reconstruction. The 
red and blue solid lines in Fig represent the tt and 7r* 
states of Si(100)-2xl AD reconstruction while the dashed 
lines represent the corresponding tt and tt* states of the 
SD reconstruction. The former clearly shows a bandgap 
~ 0.6 eV while the latter shows a continuum of states 
in the bulk bandgap region. After benchmarking these 
properties (details will be published elsewhere), the re- 
cursive surface green's function is computed for the semi- 
infinite silicon substrate. We then use a mixed-basis 
method to transfer the Si(100) surface Green's func- 
tion computed in the EHT Slater type orbital (STO) ba- 
sis into a 6-31g(d) basis set that is then connected with 
a DFT/6-31g(d) Hamiltonian for the molecule. 

While the molecule and substrate are modeled atom- 
istically, we employ a simpler treatment of the STM tip 
using a self-energy £2(^0) 0> where do is the tip to 
molecule bond length. More sophisticated models could 
be used to describe tip-sample interactions in experi- 
ments with well characterized tip structures. Vacuum 
tunneling is described using a typical WKB factor [l(| 
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FIG. 3: Energy level alignment between Ceo and silicon. 
The left panel shows the energy levels of isolated Ceo, with 
the calculated Fermi energy of the doped silicon shown as a 
dashed line. The Ceo energy levels are shifted by 1.3 eV due to 
self- consistent charging driven by the workfunction difference 
between the molecule and the substrate, and the correspond- 
ing charge transfer from Si to Ceo- Chemisorption creates 
new levels due to SiC bonds, easily seen by adding a sili- 
con dimer to the molecule and passivating the cluster at the 
bottom. Finally, including the silicon self-energy amounts to 
adding the entire silicon substrate and generates a continuum 
molecular band. The right panel shows the density of states 
corresponding to the geometry in the middle of Fig. [I] 



premultiplying £2(^0), making the net self-energy £2 en- 
ergy and distance dependant. Our approach includes the 
bias-dependent barrier profile and agrees quantitatively 
with measured STS spectra on bare silicon. 

Equilibrium band diagram. We start by describing 
the silicon-buckyball bonding chemistry, band formation, 
and the corresponding band alignment due to charge 
transfer. DFT (LSDA/6-31g) gives a good description 
of the energy levels for isolated Ceo- The highest occu- 
pied molecular orbital (HOMO) is at -6.5 eV while the 
lowest unoccupied molecular orbital (LUMO) is at -4.6 
eV relative to vacuum. Once the buckyball connects to n- 
doped silicon, electrons are transferred from Si to C60 be- 
cause the Fermi level of Si is higher than the Ceo LUMO. 
Self-consistent calculations with a Hubbard-type capaci- 
tive charging energy yield a net charge transfer that 
raises the energy levels of C60 by about 1.3 eV, which is 
very close to the difference in workfunction between Ceo 
and Si. The charging energy of the Hubbard hamiltonian 
[l^ is chosen to be 1.2 eV, consistent with experiments 
involving C60 on metal and with solid Ceo surface [2o| . 

Fig. El explains the energy level diagram using the ge- 
ometry in the middle of Fig.^as an example. The solid 
lines represent the energy levels of an isolated bucky- 
ball, while the dashed line represents the fermi level of 
bulk Si calculated from the experimental doping levels || . 
In addition to the charge transfer and band- alignment 
driven by electrostatics, there is also substantial transfer 




FIG. 4: Introducing tunnel barriers between the STM and the sample or the sample and the substrate (far left) de-emphasizes 
the HOMO levels relative to the LUMO ones due to their relatively larger barrier heights, accounting for measurements (left) 
with little or no signature of HOMO levels 0, • Although the unnormalized conductance shows this effect with increasing 
tip-sample separation d (right), normalizing the conductance (center) as in Fig. Q] restores these peaks. 



of spectral weight from Si to Ceo leading to the forma- 
tion of bonding-antibonding pairs. The energy levels of 
C60 bonded with a single surface silicon dimer system 
shows the effect of bonding, which leads to both level re- 
arrangement and level creation. A wave-function plot of 
those levels shows a lot of hybridization between C and 
Si. The rightmost panel shows the density of states of Ceo 
with the Si surface which provides the proper boundary 
conditions of the open system described with an energy- 
dependent self-energy. A clear peak appears between the 
HOMO and HOMO-1 levels in the STS (i.e., between 
peaks marked I and II), which we attribute to the Si-C 
bond arising from strong Si and C hybridization. 

Results: Peak positions and heights. Fig. E] shows the 
calculated G-Vs for Ceo docked on a clean silicon surface 
with different bonding geometries. The conductances are 
normalized using an avera ging procedure adopted in the 
experimental analyses [2l|, |22| . The upper set of figures 
corresponds to Ceo physisorbed on the Si (100)-2xl sur- 
face. The bottom of the buckyball is kept 2.1 A away 
from the Si surface dimer to ensure weak coupling. Under 
these circumstances, the STS probes the bare Ceo elec- 
tronic structure with the molecular levels (Fig. EJ) gener- 
ating G-V peaks marked I , II, III, and IV in Fig. [I] The 
vertical bars in Fig. □ denote the LDA/4-31g HOMO- 
1, HOMO, LUMO, and LUMO+1 levels of isolated C 60 
rigidly shifted by 1.2 eV due to charging as before, al- 
though the precise shift differs due to the different Fermi 
energies of the n and p-Si substrates. The excellent cor- 
respondence between the isolated Ceo levels and the con- 
ductance peaks thus gives us an elementary interpreta- 
tion of the STS data in the upper geometry of Fig. [I] 

The bonding geometry between Ceo and the Si(100) 



surface changes upon annealing from physisorption to 
chemisorption [2^ , corresponding to the middle and the 
bottom sets of plots in Fig.Q We consider two prominent 
chemisorption geometries based on experimental sugges- 
tions. The first consists of Ceo chemically bonded with 
a Si surface dimer that straddles diametrically opposite 
ends of a Ceo hexagon. The experimental dI/dV/(I/V)s 
measured using STS are reproduced by our density func- 
tional conductance calculation, with the only variable be- 
ing the geometry itself. The four main marked peaks still 
arise from the isolated C60 energy levels discussed above. 
In addition, our calculation reveals an extra smaller peak 
between peaks I and II, as in the experiments. This small 
peak has also been observed by ultraviolet photoemission 
spectroscopy (UPS) measurements on the C6o/Si(100)- 
2x1 system |2^. The origin of this secondary peak is 
from strong covalent Si-C bonding, seen in the correla- 
tion diagram between C60 and a surface Si dimer in Fig. El 
and also in the corresponding density of states. 

The bottom plot of Fig. shows the experimental 
and calculated conductances corresponding to a different 
chemisorption geometry realized upon annealing, consist- 
ing of a missing silicon dimer that causes the buckyball 
to drop into the empty trough 5]. It is clear from the 
bonding geometry (middle panel), that the closer prox- 
imity with the surface leads to the establishment of more 
covalent bonds, radically altering the electronic struc- 
ture of bare Ceo- The STS G-V in (c) is qualitatively 
different from the geometries in (a) and (b). There is 
one clear HOMO peak for negative substrate bias and a 
broadened LUMO level in the positive direction replac- 
ing the four original peaks. Our simulation captures the 
main features of this experiment. The original Si-C peak 



goes up to become the main negative bias peak, while 
several additional Si-C peaks are formed near the con- 
ducting LUMO levels due to the additional bondings. In 
addition, our calculation generates spurious peaks from 
the unrelaxed C6o structure adopted in our calculation 
for convenience. We believe that strong chemisorption 
would deform the Ceo near the bottom, eliminating these 
extra peaks by bonding-antibonding splitting. We leave 
a detailed study of Ceo relaxation for future work. 

In addition to the number and positions of the peaks, 
other experiments show variations in the conductance 
peak heights 0, 0, (Fig. BJ). A possible origin is the 
differing tip-sample spacings in these experiments. A 
WKB treatment of tunneling through varying vacuum 
thicknesses provides a qualitative explanation. Using the 
physisorbed geometry in Fig.EJa) as an example, we find 
that increasing the tip-sample gap from 6 A to 8 A and 
then to 1 nm progressively deemphasizes the role of the 
HOMO levels in comparison to the LUMO levels. Note 
however that such WKB factors are eliminated in the 
earlier log-normalized dI/dV/(I/V) plots (Fig.0J), as ex- 
pected |2l| but show up in the unnormalized dl/dV con- 
ductance plots (Fig. HJ. From the dl/dV vs V in Fig. 0] 
it is very clear that the thicker barrier cuts down the 
HOMO contributions exponentially to within the noise 
levels of the experiment. Electrons tunneling from the 
HOMO level in the negative bias direction encounter a 
higher tunneling barrier than electrons from the STM 
filling the LUMO higher up. The WKB approximation 
suggests that each peak height is reduced roughly in pro- 
portion to exp (—2kd) where k represents the decay con- 
stant of the corresponding level and d is the tip-sample 
separation. We therefore believe that measurements with 
no clear HOMO peaks are performed with an STM tip 
substantially removed from the molecule. Furthermore, 
the experiments in Fig. 0] have an extra tunneling barrier 
between the molecule and the substrate (a hydrogen pas- 
sivation layer in the top example and an oxide layer in 
the bottom [2^ ). This additional barrier would further 
de-emphasize the HOMO contributions, practically elim- 
inating them from the unnormalized conductance curves. 

In summary, we have used Ceo on silicon to demon- 
strate our capacity to theoretically deconstruct the role 
of contact microstructure on molecular conduction. The 
unexplained features are the peak broadenings, which 
should depend on coupling with the dimer and the molec- 
ular vibrational modes. A treatment of vibronic scatter- 
ing would indeed be worth pursuing both for its exciting 
physics and for further benchmarking between experi- 
mental and computational molecular conduction. 
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